HIGHER MODE DYNAMIC PLASTIC RESPONSE 
OF BEAMS WITH FINITE-DEFLECTIONS 



Carlos Antonio Pancada Guedes Soares 



Zy KNOX LIBRARY 
POSTGRADUATE school 
Y - CALIFORNIA 93940 




9 










I 













I 





































HIGHER MODE DYNAMIC PLASTIC RESPONSE OF BEAMS 
WITH FINITE-DEFLECTIONS 

t>y 

CARLOS ANTONIO PANCADA GUEDES SOARES 

// 

B.S. Kscola Naval, Lisboa, Portugal 
(1971) 

SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS 
FOR THE DEGREE OF MASTER OF SCIENCE IN CIVIL ENGINEERING 
AND THE DEGREE OF OCEAN ENGINEER 
at the 

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
May , 1976 



DimLCY KNOX LIERARY 
NAVAL POSTGRADUATE SCHOOL 




0T6ES VINHOdHVO 'A-> -x.. < 

nOOHOS 3J.VnGV!J010Gd T » \ « i 
Auvaan xo;«>i a : 



HIGHER MODE DYNAMIC PLASTIC RESPONSE OF BEAMS 



WITH FINITE-DEFLECTIONS 



by 

CARLOS ANTONIO PANCADA GUEDES SOARES 



Submitted to the Departments of Civil Engineering and 
of Ocean Engineering, on May 1976, in partial fulfillment 
of the requirements for the degree of Master of Science in 
Civil Engineering and the degree of Ocean Engineer 



A study into the higher mode dynamic plastic response 
of beams is reported here. Within the scope of infinitesi- 
mal deflections, exact solutions are obtained for generic 
symmetric and antisymmetric modes, using rigid-plastic 
analysis. In order to account for the influence of finite 
deflections, an approximate procedure is used. Also a nu- 
merical procedure which satisfies the field equations at 
different time instants is developed. Calculations are 
performed for the elastic perfectly plastic case, using a 
finite-element formulation along with a finite-difference 
time integration. Results of the different solutions are 
compared with experiments. 



Thesis Supervisor! Norman Jones 

Titlei Associate Professor of Ocean Engineering 



Thesis Reader! Jose M. Roesset 
Titlei Professor of Civil Engineering 



ABSTRACT 



1 



ACKNOWLEDGMENTS 



The author is indebted to his thesis supervisor, 
Professor Norman Jones, for his guidance and constant 
advice throughout this study. 

He is also indebted to Professor J.M. Roesset, Civil 
Engineering Department, who kindly made himself available 
to read the thesis. 

To Professor E.A. Witmer, Department of Aeronautics 
and Astronautics, who, besides allowing the use of the 
computer program JET 3 and the blast chamber for the ex- 
periments, also provided usefull advise at different 
ocasions. 

The experiments reported here were performed by 
F. Merlis and A. Abelow, and represent part of an experi- 
mental program developed by Professor Norman Jones. 

The author is also thankfull to Professor M.P. Cleary, 
Department of Mechanical Engineering and Dr. C.T. Chon, 
Brown University, for their interest and helpfull discus- 
sions. 

Professors J.H. Evans and R.W. Yeung, at different 
stages of the author stay at M.I.T. provided an invaluable 
advice and encouragement. 



2 



And, last but not the least, the author is deeply 
gratefull to his parents who provided the guidance and 
moral support during most of his student life and to his 
wife Concha who contributed significantly to his graduate 
studies by accepting patiently the often absent mind and 
providing a permanent encouragement and support which 
included the typing of the present work. 



3 



TABU-: OK CONTENTS 



ABSTRACT 1 

ACKNOWLEDGEMENTS 2 

TABLE OF CONTENTS 4 

LIST OF FIGURES 6 

NOTATION 8 

1. Introduction 11 

2. Infinitesimal Deflections 22 

2.1 Formulation of the Problem 22 

2.2 First, Second and Third Modes 24 

2.3 Symmetric Modes 26 

2.4 Antisymmetric Modes 29 

2.5 Equations of Motion 31 

3. Approximate Procedure for Finite Deflections 35 

3.1 Description of the Method 35 

3.2 Johanssen Yield Criteria 38 

3.2.1 Second Mode 39 

3.2.2 Third Mode 41 

3.2.3 Symmetric Modes 42 

3.2.4 Antisymmetric Modes 45 

3.3 Tresca Yield Criteria, Moderate Deflections 46 



4 



48 



3.3* 1 First Mode 

3.3.2 Second Mode 50 

3 . 3.3 Third Mode 50 

3.3.4 Symmetric Modes 51 

3 . 3.5 Antisymmetric Modes 53 

3.4 Tresca Yield Criteria, Large Deflections 54 

3.4.1 First Mode 54 

3.4.2 Second Mode 56 

3.4.3 Third Mode 56 

4. Approximate Numerical Procedure for 

Finite Deflections 57 

4.1 First Mode 58 

4.2 Second Mode 62 

4.3 Third Mode 63 

4.4 Numerical Timewise Solution 64 

5 . Elastic-Plastic Numerical Approach 69 

6 . Experimental Details 74 

7 . Discussion and Conclusions 78 

Figures 81 

References 95 

Appendix 1°5 

5 



LIST OF FIGURES 



Fi ,°:ure 

Figure 

Figure 

^if^ure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 



1. Free Body Diagram of a Beam Element Showing 
the Positive Directions of the Indicated 
Quantities. 



?,. Ri^id Plastic Behaviour. 



3. First, Second and Third Mode Shapes. 

4. Symmetric Modes : Half Span of the 1°\ 3 rc *, 

-th „th , . th , 

5,7, and 1 mode. 



5. 



Anti symmetr i c 

.,th ^ t h 0 th 
4,6, H 



Modes: Half Span 
and i mode. 



of 



the 




6, Tresca Yield Condition with the Circunscr i bin 
and Inscribing Square Interaction Curve. 



7 . Collapse Mechanisms 
dure . 

(a) First Mode 

(b) Second Mode 

(c) Third Mode 



Used in Approximate Proce 



8. First Mode Results • 



9. Second Mode Results- 
10. Third Mode Results • 



6 



Figure 


11. 


Result ol' the Numerical Approximate 
for the first. Mode . 


Procedure 


Figure 


12. 


Result of the Numerical Approximate 
for the Second lviode. 


Procedure 


Figu re 


13- 


Deflected Shape of Specimens 9 and 
ference [ 52 ]]. 


10 of Re- 


Figure 


14. 


Time History of the Elastic-Plastic 
of Specimen 9 of Reference [52]. 


Def lection 


Figure 


15- 


Time History of the Elastic-Plastic 
of Specimen 10 of Reference [52]. 


Deflection 



7 



NOTATION 



d Length of beam section covered wi 
h n n th time step 

t time 
t, t V q /H 

w Transverse displacement 
v C i rcunferent i al displacement 
x Axial distance 
B Beam breath 

D Dissipated energy in plastic work 
Dissipation across the i x hinge 

Energy ratio 

H Beam thickness 
1 impulse 

1 S pe c i I’ i c i m pu 1 s e 

I. Beam half-span 
M Bending moment 
M d Q H^/4, fully plastic moment 

N Axial force 

N q a Q H, fully plastic axial force 
P Pressure 



explosive 
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Q Transverse shear force 
V Vertical force, equation (4.6) 
V o Initial impulsive velocity 



W 

W* 



w 



e 



Amplitude of transverse displacement 

W 

H 

W/V 0 

WH/V^ 

Weight of sheet explosive on a beam 



W e Equivalent weight of explosive 
Wj. Final deflection amplitude 

W - . Initial Velocity, equation (3.67) 

w 1 

W Maximum permanent displacement 



a 



0 



Rotation rate across the 



. th , . 

0 hinge 



0 ^ Defined in equations (3.45,48,49,57.63) 

Y Defined in equations (3.11.17.26,33.70,71,72) 

Y Defined in equation (6.6) 

€ Strain 

h Location of plastic hinges 

-L. I. 

Location of i l plastic hinge 

n * h/L 
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n # h/v l, 

* o 

0. Rotation rate across the i hinge 
1 

K Curvature change 

X Non-dimensional kinetic energy, equation ( 2 . 13 ) 
q Mass per unit area 
P Density of' material 
a Stress 

a Uniaxial yield stress 
o J 

$ Mode shape 

Shape of the i^* 1 region of a mode 

_i_ Vn 

n„ Rotation rate across the K b hinge 

o fe < ) 

( )• fc < > 
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1. INTRODUCTION 



Since the work of Lee and Symonds [l] and of Conroy [2] 
who treated the problem of dynamic loading of rigid, perfec- 
tly plastic beams, one can identify three distint approaches 
that have commonly been used to solve the problems of dyna- 
mic loading of structures in the plastic range, namely: nu- 
merical treatments, generally including elastic and plastic 
behaviour, analytical solutions which usually neglect the 
elastic strains and the development of approximate procedu- 
res that predict only the overall characteristics of the dy- 
namic response. 

The numerical approach gives the possibility of analy- 
sing the detailed response of complex structures, involving 
both elastic and plastic behaviour (material non-linearities) 
and geometrical non-linearities such as large deflections. 
Usually these procedures can incorporate also strain harde- 
ning and strain-rate effects in the material behaviour. 

The numerical schemes available are based in a special 
discretization of the structure and a timewise approximation 
of the response. Depending on the methods used for each, 
one can consider three basic types which uses (i) finite- 
-differences in space and time [3**01 (ii) finite-elements 
in space and finite-differences in time [ 5 . 6 ] and (iii) fi- 
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nite-element discretization both in space and time L?,83. 

These methods allow not only the representation of very 
complex structures but they also provide the time history of 
the transient response. They are appropriate for use in the 
final stages of the design of a structure, when there is 
enough detailed information to provide the necessary imput 
to these numerical procedures and, when it is required to 
verufy and adjust the details of the structural response. 

However, in the initial phases, since there is not 
enough information about the structure and since one is not, 
in general, interested in the detailed response, simpler me- 
thods that provide a description of the overall behaviour, 
become more attractive. 

The main simplifications involved in the analytical 
work consist in neglecting the elastic strains, by conside- 
ring the material to be rigid for stresses smaller than a 
given yield stress. This approximation is valid when the ki- 
netic energy of the disturbance is large compared with the 
energy that can be stored elastically in the structure. It 
has been shown [9|10] that the rigid perfectly plastic theo- 
ry is a reasonable first order theory as long as the ratio 
of the mentioned energies is at least greater than three. It 
was also observed £l(>3 that for energy ratios greater than 
about ten, elastic vibrations do not have much effect on 
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the results. However, the duration of dynamic loading should 
be short compared with the natural period of vibration of 
the structure. Since most of the problems of dynamic loading 
of structures in the plastic range are produced by blast or 
impact loading and are related to the desire of having the 
structure absorb as much energy by plastic work as possible, 
before fracture occurs, very often these requirements are 
met. 

This simple rigid plastic theory has been applied to 
beams [^ 11 , 12 ], frames [_ 1 3^1 • plates [ 14,153 'and shells [l63, 
and reasonable agreement with some experimental results 
[.10,17,183 was obtained. 

However, these solutions use a small displacement as- 
sumption, in which the equilibrium equations are enforced in 
the undisturbed configuration of the structure. This assump- 
tion imposes a limit on the validity of the rigid-plastic 
idealization, since related to high levels of energy one has 
large deflections. The main developments of the rigid-plastic 
theory consisted in considering these large deflections or 
geometry changes [l9-2l3, the incorporation of strain-rate 
sensitivity [ 22-253 anc * strain-hardening [26,273 in the ma- 
terial behaviour. When these effects were included a very 
good agreement with experiments were observed [20,28,293* 

The third line of thought, recognizes that even with 
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the r igid-plast i c assumption, exact solutions are very diffi- 
cult to obtain, except for very simple situations, because of 
the discontinuities involved in the material behaviour. The- 
refore, to avoid the direct solution of the problem, approxi- 
mate procedures were developed in order to bound the final 
response or to approximate the deforming shape of the struc- 
ture. These are respectively the Bounding and the Mode Ap- 
proximation Techniques. 

The basis for the bounding techniques was laid down by 
Martin [_30] who determined a lower bound for the response ti- 
me and an upper bound for the final deflection of impulsively 
loaded rigid-plastic structures undergoing infinitesimal de- 
flections. Moralles and Neville (_3l] determined a lower bound 
for the final deflection and, Wierzbicki extended Wiartin's 
upper bound to the case of finite-deflections [ 32 ]. Martin 
[33] determined bounds for the case of time-dependent plastic 
behaviour and recently Symonds and Chon [34] developed an 
upper bound on deflection, valid for finite deflections and 
taking into account a time-dependent plastic behaviour. 

The mode approximation techniques, as presented by Mar- 
tin and Symonds [35] are based on the uniqueness proof by 
Martin [36] which in turn is an extension of an earlier work 
by Tamuzh [37]* 

Mode solutions were defined as being composed of two 
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independent'. funct i ons , one of space and other of timet 

\y = a T(t) <|> (x) (1.1) 

where a is a scalar. The structure is therefore reduced 
to a one degree of freedom system since the deformation sha- 
pe remains constant while the magnitudes change with time. It 
was shown £35] that the time function, T(t), is linear, 
which implies that in a mode solution the acceleration is 
constant. The choice of the mode shape for a given problem is 
not arbitrary but instead results from the minimization of 
the initial value of an integrated quantity that measures 
the difference in kinetic energy between the actual struc- 
ture and the mode solution. This provides an upper bound 
on the error involved. This procedure, besides it's use- 
fulness in the estimation of final deformations, provides 
also a better insight for the main features of the de- 
formation process. 

The mode approach first developed for time indepen- 
dent material behavior was extended later to rate-sensiti- 
ve material £38] and to wider types of material behaviour 
1.39]. 

Ho £4o] extended the concept to "changing mode solu- 
tions" which are solutions without the space-time separation. 
This intended to account for geometry changes and for Moving 
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or stationary discontinuity interfaces. Convergence for dy- 
namically admissible solutions was proven. It was shown that 
the changing mode solution would converge to the stationary 
mode solutions of Martin and Symonds. 

While the changing modes would correspond to the last 
few phases of the actual motion of the structure, the statio- 
nary mode solution was the last phase. 

This extension allows a better accuracy in the predic- 
tions, especially for situations where earlier phases of mo- 
tion predominate but, on the other hand, it begins loosing 
the simplicity of analysis which is in fact one of the major 
attractions of the stationary modes. 

Lee [39] extended the approach, originally for impulsi- 
ve loading, to include dynamic loading problems. Also elas- 
tic, viscous, rigid perfectly plastic and rigid viscoplastic 
behaviour were considered and it was given a wider meaning 
to mode solutions than just an approximation technique. 

Lee defined "instantaneous mode response" of the formi 

u^(x,t) = a(t) u^x.t) (1.2a) 

u^x.t) = 3(t) * i (x,t) (1.2b) 

and used the principle of virtual work to define the instan- 
taneous mode solution of specific problems. 
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He defines instantaneous elastic response as the kine- 
matically admissible field in which the potential energy is 
less than in all neighboring kinematically admissible fields 
which possess the same total displacement. When the external 
loads are zero, what corresponds to the case of free-vibra- 
tions, then the motion is in a natural mode of vibration, of 
the form i 

u i (x,t) = ^(x) T(t) (1.3) 

as the classical normal modes of vibration. 

For the rigid-plastic case, the same approach was taken 
and an analogy was made with the upper bound theorem of Li- 
mit Analysis [4l], Lee concluded that the true collapse load 
of limit analysis is the primary instantaneous mode, as de- 
fined by him. 

Therefore, one could expect that the initial motion of 
a structure would follow the primary instantaneous mode res- 
ponse if the external load is arranged in such a way that it 
is on or slightly beyond the limit load. The implications of 
this fact will be discussed later in relation to Jones work 
[47,48]. 

Recently Symonds and Chon [42] extended the mode appro- 
ximation technique to accomodate finite deflections. 

When applying this technique to problems of finite de- 
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flections, because there are travelling hinges, a stationary 
mode solution is not applicable. Symonds and Chon conside- 
red then a sequence of instantaneous mode solutions at cho- 
oser instants of time. To go from one instant in time to 
the next one a numerical iterative procedure was used, so as 
to satisfy the field equations at that specific time instant. 

Basically their mode solution satisfies all the field 
equations but does not necessarily satisfy the stated ini- 
tial velocity condition. This corresponds to neglecting the 
initial phases of the motion or, the changing modes as de- 
fined by Ho. 

However, good results were obtained and it was obser- 
ved that the shape of the structure in sucessive mode forms 
changed slowly. 

The authors noted that the numerical work in the appli- 
cation of the mode approximation, when finite-deflections 
are considered is much greater than in the infinitesimal ca- 
se. It is however substantially less than for the exact so- 
lution. Again, when extending the technique, one faces the 
problem of losing the simplicity and consequently it becomes 
less attractive. 

Recently, the bounding techniques which have been ba- 
sed mainly in energy considerations and the mode approxima- 
tion techniques which satisfy equilibrium equations after 
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having a given deformation mode choosen, seem to be mer- 
ging together. 

One of the difficulties in the mode technique consist 
in the determination of a kinematically admissible displa- 
cement or velocity field. 

Various procedures have been developed for a systema- 
tic determination of the approximate mode shape [35. etc] 
and, recently it has been shown that extremum and variatio- 
nal principles can be used to determine mode shapes [43,44], 

Independently, Jones extended the work of Sawcsuk [45, 
46] on finite deflections of plates loaded statically to 
the case of dynamic loading, developing an approximate pro- 
cedure for beams and arbitrarily shaped non-axisymmetric 
plates [47] which was later extended to shells [48]. Jones' 
method consists on the application of the principle of vir- 
tual velocities to a kinematically admissible collapse me- 
chanism. 

The main assumption is on the choice of the applicable 
collapse mechanism or displacement field. Jones assumed 
that the shape of the displacement field due to dynamic 
loads is the same as the one developed for the corresponding 
static collapse load. It was implicit on the assumption that 
the profile would be the same throughout the motion or, in 
other words, the hinges would remain stationary, although 
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plastic zones are allowed to be time-dependent. The final de- 
flection was then obtained by equating the rate of change of 
the initial energy to the rate of dissipation of energy by 
plastic work in the collapse mechanism. 

Wide comparisons with experimental results have been 
carried out [[28,29147-51] and good agreement was obtained. 

If one looks at this method in light of the so called 
mode approximation techniques, one realizes that it assumes 
a stationary mode throughout the motion. In a way, this can 
be considered a direct extension of the simple mode soluti- 
ons of Martin and Symonds [_35] to the case of finite de- 
flections. It should be noted that Jones' method keeps the 
simplicity of treatment while Symonds and Chon's extension 
[[42] is expected to give better results in the cases of si- 
gnificant changes of the hinge locations. 

The meaning of choosing the static collapse mechanism 
as the kinematically admissible displacement field can be 
established by refering to the work of Lee [_39]. This would 
correspond to choosing the primary instantaneous mode as 
being a stationary mode throughout the motion. 

Recently Jones and Wierzbicki [.52] studied theoretical- 
ly and experimentally the higher modal dynamic plastic res- 
ponse of beams, associated with the determination of defor- 
mation mechanisms capable of absorbing more energy than the 



20 



first mode. They obtained an exact solution for the first 
three modes undergoing infinitesimal deflections and used 
Jones’ method [ 47 ] to determine the influence of finite-de- 
flections on the first modal response. 

The present work consists of an extension of their 
analysis. 

Within the scope of infinitesimal deflections, exact 
solutions are obtained for generic symmetric and antisym- 
metric modes. In order to account for the influence of fi- 
nite deflections and to assess the relative merits of dif- 
ferent approximate procedures, Jones' method [47] and Sy- 
monds and Chon's method [42] were applied to the first three 
modes. Extensions to generic symmetric and antisymmetric 
modes are incorporated. 

Finally a comparison is done with a more compreensive 
elastic-plastic numerical treatment of Wu and Witmer [6] 
and with experimental results. 
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2. INFINITESIMAL DEFLECTIONS 



2.1 Formulation of the problem 

The present formulation basically follows the one of 
Jones and Wierzbicki [52]. 

A fully clamped beam of lenght 2L, subjected to diffe- 
rent forms of impulsive loading, is analysed herein. It is 
assumed that transverse displacements are small so that 
strains remain infinitesimal. 

Using the Bernoulli-Euler assumption of plane cross- 
-sections remaining plane, the displacement field of the 
beam may be approximated by the middle plane displacements. 
Moreover, if transverse shear deformations and rotary iner- 
tia are neglected, the equilibrium equations for a beam sub- 
jected to an impulsive transverse loading become » 



where the positive sense associated with these quantities is 
shown in Figure 1, and the symbols are defined in Notation. 

Combining both equations, one obtains the equation of 
motioni 



Q' - |i» = 0 



(2.1a) 



M' + Q = 0 



(2.1b) 




( 2 . 2 ) 
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The material is assumed rigid, perfectly plastic, neglecting 
therefore the effects of strain hardening and strain-rate 
sensitivity. The constitutive relations are then very sim- 
ple, as in Figure 2. For stresses smaller than the uniaxial 
yield stress, a , the strains are zero; at stress levels of 
d Q the strains become indefinitely large, of the same sign 
as the stresses; stresses higher than the yield stress are 
not allowed; 

€= 0 if I a I < a 

o 

sign €=sign a if la I = a Q 

o|s a o 

the corresponding moment-curvature relations exibit the same 
behaviour [52]. 

To complete the formulation, boundary conditions are 
imposed as appropriate for each case. 

The procedure followed is to postulate a kinematically 
admissible displacement field, introduce it in the equilibrium 
equations, use the constitutive relations and appropriate 
boundary conditions to obtain the solution. Then verify whe- 
ther the solution is statically admissible. A solution that 
satisfies all these requirements is an "exact" solution. 

In all cases considered here the velocity field is of 
the form 



(2.3a) 
(2.3b) 
(2. 3c) 
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M x, t) = 4> (x) &( t) 



(2.4) 



whore 4> if? an x-dependent function describing the shape of 



30ciated with the motion. 

Substituting (2.4) in the equations of motion (2.2.), 
results ini 



It should be noticed now that W is independent of x 
and <t> is independent of time. Therefore one can solve first 
the x-dependent part to obtain the shape and later the time- 
-dependent to obtain the motion. 

2.2 First, Second and Third Modes 

These three modes have been analysed by Jones and Wierz- 
bicki L523. They used the time-independent acceleration fields 
shown in Figure 3* These acceleration fields were introduced 
in the equation of motion, (2.2) which was then integrated 
twice spacially to obtain the moment distribution. Use of the 
appropriate boundary conditions allowed the determination of 
the acceleration and location of the hinges. 

The acceleration for the first mode was found to be 



the mode and W is a time dependent velocity amplitude an- 



M" = -p W 4> 



(2.5) 




( 2 . 6 ) 
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For the second mode it resulted 



W = 



6_ 

f in 




(2.7) 



where 



n * 2-/T 



( 2 . 8 ) 



and n L is the location of the hinge. 

It was shown that in the third mode the acceleration at 
the two hinges is the same and is given by (2.7) » where now 

n = /z - 1 (2.9) 



The permanent transverse displacements were obtained inte- 
grating equations (2.6) and (2.7) in time. The results for 
the first, second and third modes are respectively! 



W 



m _ X 
H 12 



( 2 . 10 ) 



W 



a = (/I-i ) 2 ^ 

H 6 



( 2 . 11 ) 



_b = ( y ~ 2 -\) 2 — 

H 12 



( 2 . 12 ) 



where 



v 2 2 .v 2 : 

X = ti V o L = p V o ( 2L) 

M H OH 2 

O O 



(2.13) 



is a dimensionless initial kinetic energy and V q is the peak 
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initial velocity at the interior hinge locations. 



2. 3 Symmetric Modes 



Symmetric modes are defined here as the ones in which 
the mode shape satisfies the following » 

♦ (x) = <fr(2L-x) (2.14) 

From an analysis of the first and third mode results 
£ 52 " 1 , it is apparent that only two types of shape functions 
exist in symmetric modes. 

One is applicable in the first region and is given byi 



Region li 


II 

hT 


, 0 s x s n 


( 2 . 15 ) 


corresponding 


boundary 


conditions arei 




at 


x = 0 


M = M 

0 


( 2. 16a) 


at 


x = n 
x 1 


M = - M 

0 


(2.16b) 


at 


X 

II 


Q = 0 


(2.16c) 


Substituting 


(2.15) in 


(2.5) and using (2.16), 


always 



results in one equation of motion. 

The other type of shape function isi 



W, , *1. -x 

Region ii = — 7 x — 



x-n 



i-l 



i_l sx25 i (2.17) 



w. n. -n n . -n. 

1 1 1-1 1 1-1 



where the i* h region has the i th and i th -l hinges loca- 
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ted at it's extreme points n ^ and h , and has the as- 

• • 

sociated velocity amplitudes VL , ViL_^ , as indicated in 
Figure 4. 

The applicable boundary conditions arei 



at 


X 

11 

H* 

1 

M 


M = ± M 

0 


(2.18a) 


at 


x = n . 


c 

* + 
li 


(2. 18b) 


at 


x = n i-l 


It 

O 


(2.18c) 


at 


x = n . 
1 


Q = 0 


( 2. I8d) 


these 


conditions 


in conjunction with 


equations 



(2.5) and (2.1?) results in one equation of motion for the 
i^* 1 region and in the location of the (i-l) th hinge. It 
should be noted that in the last region, one hast 

^ = L (2.19) 

and therefore these are enough conditions to determine the 
location of all the hinges and to obtain one equation of mo- 
tion for each region. 

It should be noted that equations (2.15,2.16) reduce to 
the first mode expressions [52], when n^=L , and therefore 
we can obtain immediatly the solution for the first region 
from equation (2.6)t 

W = - 6 M o (2.20) 

y n 2 

1 
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"t h 

The solution for the i region, is obtained by subs- 
tituting (2.37) in (2.4). 

IJso of (2.1B) leads to 



W. 

i 




w 



(n . -n 

' i 




12 

4 




( 2 . 21 ) 



( 2 . 22 ) 



• # 

where W is given by (2.20). 

Introducing (2.20) in (2.22) gives the general expres- 
sion for the location of the l hinge » 



* i = n i _ 1 + ^2 n 1 , i>2 ( 2 . 23 ) 

Now, when we are considering the last region of each 
mode shape, equations (2.19) and ( 2 . 23 ) result in 

L = n + ’ / T h 1 (2.24a) 

and also from ( 2 . 23 ) 

^ i _ 1 = B i-2 + *^2 n i (2.24b) 

Equations (2.24) lead toi 

^2 = L -2 * r 2 h 1 (2.25) 

sucessive substitution of ( 2 . 23 ) in ( 2 . 25 ) leads to the ge- 
neral form* 

n Q = 0.0 (2.26a) 
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n = b-i ^2 n i = 0,1,2, .. 

n-i ± t 


. ,n (2. 2ob) 


where n is the number of interior hinges per 
the (2n-l) th mode. 

When i = n-1 , (2.26b) becomes i 


half span of 


n = L 

1 l+(n-l) ^2 


(2.27) 



Substitution of (2.27) in (2.20) and (2.21) gives the 



general expression for the acceleration, valid 
ons « 


in all regi- 


W = - 6 5[l+(n-l)^ 

4 L 

and (2.26b) becomes! 


(2.28) 


Li = H-li-l fz , i - x ,2 n 

L l+(n-l)/T 


(2.29a) 


o 

II 

o 


(2.29b) 


It should be noticed that for n=l and 2, 


equation 



(2.29) gives the previously obtained (. 52 ] hinge location for 
the first and third mode 

2.4 Antisymmetric Modes 

The antisymmetric mode shapes satisfy the following re- 
lation! 
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♦ (x) = - $(2L-x) 



( 2 . 30 ) 



In these modes there exists three types of shape func- 
tions ♦j. Two of them are given by equations (2.15) and 
(2.17) of the symmetric modes and, the third type corres- 
ponds to the region just adjacent to the midspan of the 
beams 

Vl “Hr • V XSL ( 2 - 31 ) 

n 

the applicable boundary conditions ares 



at 


tr c 

II 

X 


• 


M = ± 

0 


( 2 . 32 a) 


at 


x = L 


» 


M = 0 


( 2 . 32 b) 


at 


X 

II 

3 


f 


Q = 0 


( 2 . 32 c) 



where the second condition results from antisymmetric con 
siderations. 

Introduction of (2.31) in (2.5) and satisfaction of 
(2.32) results ins 




Introducing equation (2.23) for i = n in equation 
(2.33) results ins 

Vi = L - (| + !) /Thi (2.34) 

sucessive substitution of (2.23) in (2.34) results ini 
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(2.35) 



n . 

n-i 



L - 




i) 




1 



when i = n-I this expression will give the location of 
the first hinge after the supports, as a function of mode 
number i 




1 

l+» r 2 (n- A) 



( 2 . 36 ) 



Introducing this in equations (2.20) and (2.35) will 
result ini 

M 

w = - £— ! [1 ♦ (n- h «^2] 2 (2.37) 

p L 



and 



h 

o 



0 



(2.38a) 



i = l+/z (i-1) 
l+*2 (n- \) 



i = 1,2, . . .n (2.38b) 



which are the general results applicable to the 2n x mode. 

For n = 1 equations (2.37) and (2.38) recover the 
results of the second mode, equations (2.7) and (2.8) as 
obtained in [ 52 ]. 



2. 5 Equations of Motion for Symmetric and Antisymmetric 
Modes 



As resulted from the previous analysis of symmetric and 
antisymmetric modes, the equation of motion that governs all 
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the modes is the same i 



W 




(2.39) 



where h ^ is the distance of the first hinge from the clam 



ped end of the beam, and assumes different values, as given 
by equations ( 2 . 27 ) and (2.36). 

It should be noticed that for a given mode, it results 
from (2.3a) that the acceleration is constant. Indeed, Mar- 
tin and Symonds [_35]showed that for stationary mode soluti- 
ons, like the present one, the velocity varies linearly 
with time, which is the same as having a constant accele- 
ration. 

The initial conditions applicable to the cases conside- 
red here are, at t = 0, 



what corresponds to applying an initial velocity of trian- 
gular shape and maximum value V q to an undeformed beam. 
The shape of the initial velocity imparted to the beam is 
equal to the corresponding mode shape as determined from 
the previous analysis, equations (2.15» 17* 31) and Reference 



W = 0 



(2.40a) 



W = V 



o 



(2.40b) 



[52] 
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Integrating (2.39) twice with respect to time and ap- 
ply i nr the initial conditions (2.40), results in the deflec- 
tion amplitude i 



W 



V t 
o 




(2.41) 



The final deflection occurs when the velocity vanishes. 
This will happen at* 



[i V n ? 

+ - o 1 

x f 6 M Q 

and substituting in (2.41), we will obtain the final 
tion amplitude: 



( 2.42) 



def lec- 



V _f . V o P 1 
H 12 M q 

which can be rewritten as: 



(2.43a) 





(2.43b) 



where X is given by (2.13)» 

As the mode number increases, the corresponding value 
of decreases and, consequently, the final deflections 



decrease also. 

The final deflected shape of the beam is given by: 
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(2.44) 



whore ^ 
For 
obtained 



W, 



1.2 



(r 1 ) 



*(x) 



x) in the appropriate mode shape, 
the first three modes (2.43) <tives the 
by Jones and Wierzbicki [.523. 



results 
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3. APPROXIMATE PROCEDURE FOR FINITE DEFLECTIONS 



3.1 Description of the Method 

The approximate procedure used here is based on the the- 
oretical study of Jones [4y] mentioned previously. It con- 
sists of assuming a kinematically admissible collapse me- 
chanism and then, by using the principle of virtual veloci- 
ties, equate the energy imparted to the structure to the 
energy dissipated by plastic work at the hinges. 

The procedure was developed for a general non axis- 
symmetric, initially flat, rigid perfectly plastic plate and 
a specialised form was presented for a plate which deforms 
into a number of rigid regions separated by straight line 
hinges. 

The balance of energy for the particular case discussed 
above can be expressed by [47,52]» 

S (P-|iw) wdA = E S (Nw-M) 0. dZ (3.1) 

A i=l Z 1 

where the area A on the left hand side of the equation is 

the total area of the plate which is subjected to the exter- 

. • 

nal pressure P, 0^ is the relative angular rotation rate 
across the i u hinge and Z is the length of the hinge 
which, in the case of the beam is it's width. 
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The energy dissipated per unit width of the beam is gi- 



ven byi 



n 



D = E 
i=l 



(Nw-M) 0 i 



(3-2) 



while the left hand side of (3-1) represents the work done by 
the external and inertia forces. 

Equation (3*1) can be specialized for the case of a beam 
under impulsive loading by setting P=0 and performing one 
integration across the width of the beam, resulting ini 



In the present study we will therefore apply equation 
(3.3) to the different modes. 

The imediate consequence of dealing with finite-deflec- 
tions for structures axially restrained consists on the de- 
velopment of membrane forces which increase in magnitude as 
the deflection increases, and provides a considerable reserve 
of strength. 

In passing, it should be noted that great care should 
be taken when modeling structures as axially restrained sin- 
ce Jones showed for the static case [53] that remarkbly small 
in-plane displacements at the supports can change the respon- 
se from that of an axially restrained beam, with considerable 



S - p w w dx = E (Nw-M) 0 . 
x i=l 1 



i=l 



(3-3) 
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reserve « trengthen i rig beyond the limit load for moderate la- 
teral deflections, to that of a freely supported beam with no 
increase in strength beyond the limit load. 

Neglecting the effect of transverse shear forces L5*0as 
is appropriate at least for the lower modes [_52], we must 
deal with the interaction curve that relates the moment and 
axial force. 

The Johanssen yield criteria, sometimes known as the 
square yield curve, is known to represent well to behaviour 
of concrete structures 5571 • 

The Tresca and the Mises-Hencky yield criteria repre- 
sent better the behaviour of metal structures, and while the 
second is more accurate and more suited for numerical work 
[6], the Tresca yield condition being piecewise linear is 
in general choosen for analytical work. 

However, Johanssen's yield condition is still much sim- 
pler to use than Tresca and, is has been shown (_ 25, *+7,563 
that even for metal structures, a square yield curve cir- 
cunscribing Tresca' s curve could be considered to provide a 
lower bound to the maximum permanent transverse displacements, 
while an inscribing yield curve, as shown in Figure 8, (with 
sides 0.618 times smaller) could be considered to give an 
upper bound. 

Therefore, in the present study both yield conditions 
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will be used and the results compared. 

In all cases, the assumed collapse mechanism has the 
same shape as the corresponding infinitesimal deflection re- 
sult, and because of symmetry only half-span is analysed. 

3.2 Johanssen Yield Criteria 

Using this criteria, the yield condition and normality 

requirements will be satisfied when one has at the hinges 

N = N and M = ± M . 
o o 

It should be noted that for moderate deflections, in 
the absence of longitudinal accelarations , equilibrium re- 
quires the axial force to be constant throughout the beam 
[ 5?] and, consequently in the present case equal to N . 

In solving finite-deflection problems, very often a cons- 
tant axial force is assumed [_ 19 » 25 »etc] . 

It can be shown that the left hand side of equation 
(3*3) for half a beam is given by« 

L 

f -pwwdx=- — WW (3» i +) 

o 3 

for all modes. 

Now, for each mode, the rotation rates should be deter- 
mined and introduced in (3*2) to determine the dissipated 
energy which, when equated to (3.4) will result in the 
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equations of motion. 

This procedure has been used by Jones and Wierzbicki, 

[<S?] who applied it to the first mode. They obtained for the 
maximum deflection! 

-~ 2 - = — -{(1+X/3) 1/2 - 1} (3.5) 

Here we will apply it to the modes higher than the first. 

3.2.1 Second Mode 

The velocity field is given by equation (2.4) where the 
mode shape, shown in figure 7b , is the same as for infini- 
tesimal deflections [52] 1 



region 1 1 


$-1** * 

y l n 


VI 

X 

VI 

0 


( 3 . 8a ) 


region 2 1 


L-x 

*2' L-r 


n s x s i, 


(3:8b) 



The corresponding rotation rates are 1 



at x = 0 


1 

d = 


. Jl 


( 3. 7a) 






n 

• 




at x = n 


a = 


W L 


(3.7b) 



n (L-n) 

where q is given by equation (2.8), from the infinitesimal 
deflection analysis [52]. 

The normality requirements are satisfied when 



at x = 0 


N = N , 


M = M 

0 


(3.8a) 


at x = q 


N * N , 

0 ’ 


M =-M 

0 


(3.8b) 
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If (3.7) and (3.8) are introduced in (3*3) we will 



obtain 1 



where 



and 



W* + y 2 W* + = 0 



<• 




W H 





W, 



W 

H 






l!o 

H 




12 12 

X * 3 n/2-4 



The initial conditions are, at t* = Oi 



(3-9) 

(3.10a) 

(3.10b) 

(3.10c) 
(3. 10d) 

(3.11) 



W* = 0 (3.12a) 

w* = 1 (3.12b) 

Solution of (3*9) and use of (3*12), in conjunction with 
the fact that the final deflection occurs when W* = 0 gives 
the duration of motioni 

t* f = - tan -1 { 2 -^) (3.13) 

Y Y 
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and the final deflectiont 



W» f = — ” { (1+X - ^ - 1) (3.14) 

2 \F? 3 

Results are shown in Figure 10. When X (6 \[ 2 - 8)/3 «1 
equation (3*14) reduces to the corresponding infinitesimal 
case result [52], equation (2.11) here. 

3.2.2 Third Mode 

As shown in Figure 9c, the velocity field is given by 
(2.4) where [52]» 

region 1* <J>^ = ^ , 0 < x s' 1 (3.14a) 

w 

region 2i i ^ , hs x s L (3.14b) 

^ # 2 L-n L-n 

The rotation rates (3*7) remain valid in the present 
case and, in addition we have t 

at x = L n = - (3.15a) 

where now, from the infinitesimal deflection [52], n is gi- 
ven by (2.9). 

Also, in addition to the normality requirements (3.8) 
we have i 

at x = L N = N o , M = M 0 (3.15b) 
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the resultant equation of motion is given by i 



W* + + 7 X j= = 0 

* * 6 - 2/2 

where ( 3 . 10 ) remains valid and now, 

2 _ 12 . 2 - 3/2 
Y X 4-3/2 

This equation predicts the final deflect i on 1 

w,f = -ip {(1+ k S0^P ) 1/2 - x, 

6 - 2/2 3 2 - 3/2 



(3. 16) 



0 . 1 ?) 



( 3 - 18 ) 



which occurs at 




( 6-2/2 j 



(3.19) 



Results are shown in Figure 10. When X(80-57 '/2)/(6-9 sf 2) «1 
equation ( 3 * 19 ) reduces to the corresponding infinitesimal 
case result [ 52 ], equation ( 2 . 12 ) here. 



3.2.3 Symmetric Modes 



The collapse mechanism has the same shape as given in 
section 2.3* The axial force is N q throughout and normali- 
ty is satisfied having at the support M = M q and in conse- 
cutive hinges M = + M q . 

Equation (3*4) remains valid. To complete the formula- 
tion, we must now determine the dissipated energy, for a 
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J. u 

generic 1 hinge, an shown in Figure 4. 

The rotation rate across the hinge, is g i von byi 



9. = 2W 



<Vi - n i-i > 



<r 'i- r * l -i ) (n iti- n i ) 



(3.20a) 



t hi 

The location of the i hinge in the symmetric modes, is 
given by (2.29) and, if this is substituted in (3.2.0a) we 
obtain i 



0. = [1 + \[ 2 ( n-1) ] (3.20b) 

1 L 

The dissipation function corresponding to the i^* 1 hinge is 
then given by« 

v/2 M n W w 

D = 2 — (1+4 -) [l+>/F(n-l)l (3-21) 

1 L H 



where use was made of equation (3*2). 

In the first region, which includes the hinge at the 
support and the first interior hinge, the dissipation func- 
tion is given byi 

D x = M q ( 1 +(W2) (1+4 g)> (3.22a) 

or, using (2.29), 

NI W ... 

D 1 = -~—(l+(l +>f2)(l+4 §)>Ll+(n-l)/2] (3.22b) 

At the n^ hinge, which is located at x=L, the dissipa- 
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tion function is given by » 



D = ^ D. 
n 2 i 

since, only ho If of the rotation is included because 
half of the beam is considered. 

"t h 

The total dissipation for the (2n-l) T,! mode is t 



by » 



D = D 1 + (n-3/2)D i 



or 



MW ■) 

D = — — fl+(l+4 -) [l+2/2(n-l)] [l+(n-l)v/2] 
L H J 



Equating: (3.24) to (3*4) results in the equation 



tion t 



vJ + J- w + x! = 0 

* Y * 2 1+2/7 (n-1) 



where 

Y 2 = ~ (_ 1+2/2 (n-1)] |_l + (n-l)/2] 
The duration of motion is given by 
t. f = ^ tan' 1 
and the final deflection isi 



W 



3 I 1+J2" (n-l)P 



*f " 2+ 



L1+/F (n-1)]- 



(3.23) 

only 

i ven 
(3.24a) 

(3.24b) 
of mo- 

(3.25) 

( 3 . 26 ) 

(3.27) 

( 3 . 28 ) 
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For n=l and n=2 this result reduces to equations 
(3.5) and (3*18) of the first [ 52 ], and third mode respec- 
tively. 

3.2.4 Antisymmetric Modes 

As can be observed in Figure 5» .in these modes only 
two types of dissipation function exist. One is valid for 
all but the region closer to the support and the other type 
is valid for the first region, containing the hinge at 
the support and one interior hinge. It is given byi 

mw r ~ 

D x = -2^- [ 1+ ( 1+4 |[) ^|^][l+/2(n- |)] (3.29) 

if use is made of equation ( 2 . 36 ). 

The rotation rate across the i u hinge is given by 
(3* 20a). Using this equation and ( 2 . 36 ), one can derive the 

4* Vi 

applicable dissipation function for the 1 hinge 1 

2j2 M n VI w _ , 

D i = xT — ( 1+4 t 1+ & (n- f)] (3.30) 

4 * 

The total dissipation for the 2n mode is given byi 

D = D 1 + (n-1) D i (3.31a) 

or i 

M W 

D = ( l+( 1+4 ^)[l+2v/2(n- -|)] } [l+y/2(n- |)] (3. 31b) 
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Equating (3*31) to (3.^) results in the equation of mo- 



tion i 



W + W + . = o 

* Y * 2 l+2\f2(n-l/2) 



(3.32) 



where 



= ~ [1+/2 (n-1/2)] [l+2^(n-l/2)] 



(3.33) 



the final deflection is given byi 



W *f = (t 1+ 3 



lH-2/2( n-1/2) -j 
[l+j2(n-l/2)] 3 




and it occurs at 



'*f 



= — tan 



-1 



$ 






(3.35) 



When n = 1 these equations reduce to the second mode re 
uults, equations (3*9) to (3.13). 



3.3 Tresca Yield Criteria - Moderate Deflections 

In accordance with this criteria, the interaction curve 
between axial force and moment is parabolic, as shown in 
Figure 8 and can be expressed by \_9~]x 

2 

+ -\ = 1 (3.36) 

N o 

The normality of the strain-rate vector to the yield 
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surface in satisfied by the associate flow rulet 



M K N 

-2-r = — for |N| < N (3.37a) 

N € 2N 

o 



, M ic , 

•^S — — r for |M|= N (3.37b) 

^ N C 

o 

• • 

where € and K are the generalized plastic strain rates, 
the extensions] strain rate and the curvature rate respec- 
tively. Using this criteria, Jones [47] derived the dissipa- 
tion function for the interior hinge of a simply supported 
beam deforming in the first mode » 



D. 

1 



M (1+4 9. 

o h 1 



(3-38) 



For the case of a clamped beam, it becomes 1 

2 . 

D. = m ( 1+3 e. (3-39) 

X U 1 

Also in a clamped beam, one has for the hinge at the support! 

D i = M o (1 - «i <3- 40 > 

rl 

These dissipation functions are valid for N<N . 

o 

When N = N q , it can be seen from equation (3*36) that 
the moment must vanish. When that happens, these equations 
are no longer applicable. It can be shown that for a clamped 
beam that occurs when W # = 1. 
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The present analysis are therefore valid for deflections 
smaller than the beam thickness which are called here modera- 
te deflections. 

In all the cases only half of the beam is analysed be- 
cause conditions of symmetry or antisymmetry apply. 

3.3.1 First Mode 



The velocity field is given by (2.4) where the shape is 
the same as used by Jones and Wierzbicki [_52] for the infi- 
nitesimal case and when using the previous yield criteria. 

It is shown in Figure 8a and is given by» 



<p = ■£ , OsxsL (3.410 

the corresponding rotation rates are t 

at x=0 9 = - | (3.42a) 

at x= L a = ~ (3.42b) 



Using the appropriate rotation rates in conjunction 
with (3*39) and (3.40) and substituting all in (3*3) » re- 
sults in the equation of motiont 



W + 



6 M 

o 





6 M 
o 




(3.43) 



which agrees with Jones [^4?]* when his equations of motion 
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are specialized for the case of impulsive loading. 

This equation can be rewritten in the non-dimensional 

formi 

W* + W 2 = e 2 (3-44) 



where (3* 8) holds andi 



3 




3 



2 




(3.45a) 

(3.45b) 



the initial conditions to be satisfied are given by (1.12). 



Solving this equation by sucessive approximations, one 
obtains when using the second approximation and condition 
( 3 - 12 ). 



w * = t, + 



S 2 



t ^ 



12 



33 ? t* 3? t* 

(1+ — | — + -~5 — ) (3.46) 



and 



= 



Pi 33 ? t* 33| t| 

1+ »2 ** " “TB~- (1+ -T“ + “To— } 



(3.47) 



the duration of motion is obtained from equation (3*47) when 
satisfying the condition W # = 0. Substitution of the smal- 
lest root of the resulting equation in (3*46) will give the 
final deflection amplitude or, the maximum deflection. 
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3.3.2 Second Mode 



It can bo shown that for this collapse mode, described 
in 3.2.1, the dissipation function applicable to the interior- 
hinges is still given by ( 3 • 39 ) • Using this equation and 
( 3.40) , the equation of motion obtained is still (3*44) whe- 
re now 1 



a . 1 JtzM 

1 x 3 ,/?-<» 



(3.48a) 



»2 




2 

3/2-4 



(3.48b) 



the solution for this mode is obtained from (3.46) and 

(3.47), using the new values of 0. 

1 • 



3.3.3 Third Mode 

The same velocity field of section 3*2.2 applies here. 
The dissipation function for the hinge at the support and 
the hinge next to it is, as in the second mode, given by 
equations (3*39) and (3*40). For the central hinge, equation 
(3.38) is applicable. 

The equation of motion is given by (3*44) withi 



_ 6 6/2-6 

1 X 3/2-4 



(3.49a) 
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( 3,49b) 




6 >g; 

x 3/2-4 



the final deflection is obtained from (3*46) and (3*4?) • Re- 
sults from ail these modes are shown in Figures 8-10 com- 
pared with the solutions obtained using the previous yield 
criteria. 

3.3*4 Symmetric Modes 

The total dissipation in these modes is obtained by 
considering all the three types of disspation function 
equations (3.38-3*40), It is given byi 



For the first mode (n=l) since the central hinge is at 
the same time the first after the support, we havei 



D q is the dissipation function at the support and 
is given byi 



D = D + D, + (n-2) D. + D 
o 1 1 n 



n ^ 2 ( 3 . 30 a) 



for the (2n-l)^ mode. 



D = D + tt 
o 2 



0.50b) 



D o = a- 4 ) Z U+/2(n-l)] 
n 



( 3 . 51 ) 



where use was made of (3.40) and ( 2 . 27 ). 
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A. U. 

For the i hinge the rotation rate is given by (3.20b). 
Using this in conjunction with (3*38) and (2.27) results in 

2 * 

D. = (1+4 %) - 2</2 |_lW2*(n-l)] (3»53) 

1 Yi L 

The n th hinge occurs at mid-span and the applicable dissi- 
pation for all modes except the first is given by 



D 



n 




(3.54) 



Therefore, for all but the first mode, combining (3*50) to 
(3*54), one obtainsi 



D = 



M W 
o 



[2+>f2(8n-9)]+[2+2^f2(n-l)]} j_l+x/2(n-l ) ] , n*2 (3-55) 



Equating (3.55) to (3*4) results in 

w * + wj - 0 2 (3,56) 

where 

0 1 = X L2+/2(8n-9)] Ll + ^2(n-l)] (3. 57a) 

6 r- 2 

0 2 = " X [l+V2(n-l)] (3.57b) 

The solution to this equation is obtained from equations 
(3*46) and (3.47) using the present values of 0 . When 
n = 3 equations (3-57) reduce to (3.49). 
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3.3*5 Antisymmetric Modes 



The dissipation function for the 2n^ h mode is given by» 



D = D o + D 1 + (n-1) D i (3*58) 

where now 

= (1- ~) - (_l+^2(n-l/2)] (3*59) 

0 n L 

D, = (1+3 ^ ) - (1+/?) Ll+/2(n-l/2)] (3*60) 

1 L 

D. = (1+4 - 2 \fz [l+/2(n-l/2)] (3*61) 

1 L 



Combining these equations will obtaim 
MW 2 

Da- 2 — {(_2+^2(8n-5)] ~~2 ♦ L2+\f2( 2n-l)] J[l+/2(n-l/2)] (3*62) 

L H 

Substituting (3.82) in (3.3) and using (3*4), will re- 
sult in the equation of motion, equation (3*58), where now 

* 1 = £ [2+/2(8n-5)] Ll+\f2(n-l/2)] (3.83a) 

0 2 [2+j2(2n-l)] |_l+vf2(n-l/2)] (3.63b) 

the solution of the equation of motion is given by (3*46) and 
(3*47). 
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3.4 Tresca Yield Criteria - Large Deflections 



The results obtained in the previous section are valid 
for deflections smaller than one beam thickness. For greater 
deflections, the beam exibits large deflections, as defined 
here, in which case it behaves like a plastic string. This 
has been reported elsewhere [|e.g. 19,47]. 

We will analyse here only the first three modes since 
it is not likely to have modes higher than these enter this 
range of deflections. 

3.4.1 First Mode 



The 


dissipation function applicable in this case, as 


given by 


Jones [.47] is» 




D i = 4 M q H f>i , H » 1 (3.64) 



Equating this to (3*4), taking (3.42) into account, one 
obtains the equation of motioni 





W* + y 2 W* = 0 (3.65) 


where 


Y 2 = p (3.66) 


The 


initial conditions for this phase of motion are gi- 
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1. 



von by equations (3*3*0 to (3*36), when = 

The duration of first phase is given by the smallest 
root of the equation resulting from substituting W # = 1 in 
equation (3.46). Substituting the resulting value of time 
in (3.47) will result in the initial velocity for this phase 
of motion, W^. . 

Then the initial conditions are at t^ = t # . 1 

W* = 1 (3.67a) 



W„ = W„. 



(3.67b) 

Solving equation (3*65) subjected to these initial con- 



dition results in 1 



W^. 



w. 



* T 

= -rr siny(t*-t #i ) + cosy(t # -t #i ) 



(3-68) 



the final deflection occurs when the velocity vanishes ati 

W„ . 



'*f 

and the final deflection isi 

.2 

W*. 

w* f = (1+ -r) 

Y 



“1 -M- ^ 

= t* ? + tan ( —■ ) 

0 



(3.69) 



1/2 



(3.70) 
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3.4.2 Second Mode 



Here, equations ( 3 * 65 ) and ( 3 . 68 ) to (3.70) of 
previous analysis, are still valid. However, equation 
(3.66) becomes now 

2 = 12 12 
X 3^2-4 

3.4.3 Third Mode 

The analysis remains the same as for the previous 
except that now the different rotation rates, lead to 

y 2 = T ( 5 + 3V2) 



the 



(3.71) 



modes 



(3.72) 
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A„ APl’HOX 1 MATE NUMERICAL PROCEDURE FOR FINITE DEFLECT IONS 



In the cane of infinitesimal deflections, considered in 
Reference [_52] and i- n section 2 we were able to obtain an 
exact solution with a stationary mode shape because the 
initial velocity was arranged so as to have exactly the same 
shape as the mode solution and, because only bending moment 
was considered, consistent with the infinitesimal range of 
deflections. Only when these two conditions are satisfied 
simultaneously can an exact solution with a stationary mode 
be obtained. 

For example, changing slightly the loading systems to 
be a uniform initial velocity, instead of triangular, will 
result in having one phase of motion in which two hinges 
leave the supports and travel inwards. Only when the two 
hinges meet together will the mode shape considered in Refe- 
rence [ 52 ] or equation (3*^1) prevail. 

Therefore a solution only with this shape as being 
stationary would only approximate the motion and therefore 
would not be exact. 

The inclusion of finite deflections in analysing cases 
where there are axjal restrains, changes also the nature of 
the problem and stationary shapes can no longer be considered. 
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Therefore hero wo will consider velocity fields of l.lm 



formi 

w (x,t) = W( t) <J> (x ,t) 


(4.1) 



where the shape of the mode is allowed to change with time. 
The equilibrium equations for finite deflections are 



[57]. 




»• » 1 

M - ( N w ) + p w = 0 


(4.2a) 


Q = - M* 


(4. 2b) 


n’ = 0 


(4.2c) 


in the absence of transverse and 


longitudinal pressure loa- 


ding. From (4.2a) and (4.2c), it 


results i 


" " 

M -Nw + p. w = 0 


(4.3) 


It is apparent from section 


3 that the consideration of 


the Johanssen or of Tresca yield 


condition do not change the 



general concepts involved and that in general Johanssen's re 
lation leads to less complicated treatments. Therefore, in 
this section, only the square interaction curve will be used 
to obtain the results. 

4.1 First Mode 

As now we will allow the mode to change with time, we 
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will not restrict the interior hinge to be located at cen- 
ter of the beam. Therefore, we will formulate the problem 
for the entire beam, of length 2L, while in the previous 
analysis only half span was considered. 

We will however be restricting the solution in that we 
will not account for situations such as two hinges spreading 
from the center of the beam, leaving a rigid region between 
or, any other mode form, except the one which has only one 
hinge in the interior of the beam. 

The mode shape is then given byi 



region 1 * 4> ^ (x ,t) = - 



0 s x ^ h 



(4.4a) 



region 2i <f> 2 (x,t) = 



n £ x £ 2L 



(4.4b) 



where now h varies with time and its time derivatives must 
be considered. It should be noted that this mode shape is 
not a symmetric one and is therefore applicable to any tri- 
angular initial velocity condition, with it’s peak at x=n. 

Equation (4.2c) together with the yield criteria impo- 
ses that N = N q throughout. 



The 


applicable 


boundary conditions are i 




at 


x = 0 


M = M 

0 


(4.5a) 


at 


x = n 


M = - M 

o 


(4.5b) 


at 


x = 2L 


M = M 





o 
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at x = n 



V = 0 



(4, 5d) 



where V is defined as the vertical force and is Riven byi 



V = Q+N w ' 
o 



(4.6) 



Introducing now (4.1) and (4.4a) in (4.3) and using 
conditions (4.5a,b,d) results ini 

3 N 6 M 

w - Aml - + — % w + — % = o (4.7) 

n n ^ (1 ^ 

Use of (4,4b) and (4.5b-d) in the equilibrium equation, 
(4.3) results ini 



.. w 4 

W - + 



3 N, 



6 M 



2L-4 q ( 2L-4 ) 2 q ( 2L-4 ) 2 



= 0 



(4.8) 



• I 

If one now introduces W from equation (4.8) in (4.7), 
the result isi 



( L-4 ) 2 6 M o 

2L- 4 qW 4 l 



( 1+2 2 ) 
H 



(4.9) 



this equation gives the velocity of propagation of the hinge. 
It should be noted that it vanishes when ri = L. 

The applied velocity field is symmetric and has it's ma- 
ximum value at x=L. Then, the initial conditions are at 
t = Oi 
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n = L 



(4.10a) 
(4. 10b) 



W = 0 

W = V (4.10c) 

o 

Therefore the hinge is initially at n= l and the velo- 
city of it's propagation clearly vanishes. In this case we 
have therefore a stationary hinge as given by this formula- 
tion. 

Therefore equations (4.7) and (4.8) reduce just to 

.. 3 N n 6 A/I 

W + |w+ = o (4.11a) 

p h |ili 

It's non-dimensional form is s imply i 

W* + ~ W* + £ = 0 (4.11b) 

It should be noted that because of the stationarity of 
the hinge, this solution reduces to the case considered by 
Jones and Wierzbicki [52] when taking finite deflections in- 
to account. 

Also, the stationarity of the hinge is the result of the 
symmetry of the initial velocity. The formulation, equations 
(4.8) and (4.9) handles the more general case of non-symme- 
tric initial conditions, in which case equation (4.9) pre- 
dicts a moving hinge. 
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4.2 Second Mode 



The velocity field is given in this case by (4.1) and 
(3.6) where in the last one n varies with time. The bounda- 
ry conditions are given byi 



at 


x = 


0 


M = M 

o 


(4.12a) 


at 


x = 


n 


M = -M 

0 


(4.12b) 


at 


X = 


L 


M = 0 


(4.12c) 


at 


X = 


n 


V = 0 


(4. 12d) 



where (4.12c) results from antisymmetry and V is given by 

(4.6). 

Introducing the velocity field corresponding to the 
first region in the equilibrium equations, (4.3) and using 
the appropriate boundary conditions results ini 

.... 3 N 6 M 

w w + % =0 (4.13) 

n t i n 2 qh 2 

Performing the same for the second region, results int 

* 3 N 3 M 

W + W — + 2-s W + = o (4.14) 

L-n (.i ( T. Jn ) 2 p(L - n ) 2 

Now we must solve these two equations for the deflec- 
tion W and for the hinge location q . In non-dimensional 
form they will be respectively* 
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(4.15) 



n 

* 



6 

xnTW 






12 W # 

rnX 



C" 



w. 



3(2-%) 12 w, 

X X %(!-%) 



(4.16) 



where 



n * 




L 



(4.1?a) 



n 

* V L 



(4.17b) 



o 

These equations will be solved numerically, as indicated la- 
ter. 



4. 3 Third Mode 

The velocity field is given by (4.1) and (3*14) where 
h is to be understood as time varying. The applicable boun- 
dary conditions are« 



at 


X 


= 0 


M = M 

0 


(4. 18a) 


at 


X 


= n 


M = -M 

0 


(4.18b) 


at 


X 


= L 


M = M 

0 


(4. 18c) 


at 


X 


= n 


V = 0 


(4. I8d) 


at 


X 


= L 


V = 0 


(4. I8e ) 
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Introduction of the velocity field (3.l4) in the 
equilibrium equations, (4.3), and use of conditions (4.lMa-c) 
leads to the expressions of moment in the two regions. Use of 
conditions (4,l8d,e) will result in the equations of motion, 
which can be written ast 



^•1 



8 

X n * (1- n*) 






X ( 1-n ) 1 



(429) 



W *o = 



_8_ 

n 



w „ +w„.n * 

*0 *1 * 



(1- 



n } 2 
* / 



- (w* 0 +w #1 ) 



+ w 



1- 



* 1 



( 4. 2o ) 



w # n # 1? w„ 

n - + SSL 



w 



*0 



* n * 



* w *0 n * 



(4.21) 



# 

where now the velocities in the central region and on the 
lateral regions W q are no longer equal. 

These equations will be solved numerically. 



4.4 Numerical Timewise Solution 

While for the infinitesimal case it was possible to 
obtain an analytical solution of the equations of motion, 
now we will solve them numerically since they don't seem 
amenable to a closed form solution. 

There exists a wide range of integration operators 
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that have been unoci sucessfully in structural dynamic pro- 
blems L5*0* Their characteristics as to stability, converge 
and damping are well known [59]» except for the recent and 
apparently more attractive Park's method L60,6l], 

Because in the present problem the elastic vibrations 
are neglected no major stability problems are expected and, 
instead of the more sofisticated operators such as Park and 
Houbolt , a Newmark 0 operator [ 62 ]will be used. 

The Newmark 0 operator, for y= 1/2 and 0 = 1/4, also 
known as the trapezoidal method or the average acceleration 
method [ 63 ] » exibits very good characteristics for linear 
problems, such as being unconditionally stable (0 = 1/4) and 
introducing no false damping (y = 1/2). 

However, for non-linear problems, stability is not 
assured. 

The average acceleration method is an implicit method 
in that it needs information about the time step to be cal- 
culated, before it calculates it. Basically it uses for the 
value of acceleration during one time step the average value 
of the accelerations at its begining and end. Therefore, if 
no extrapolation scheme is used, as was choosen here, itera- 
tion becomes necessary. 

In general, for the n l time step, the values of the 



65 



velocities and deflections are given by« 






(W +W . ) 
n n-1' 



(4.22a) 



W 



n 




(w +w 

n 



n-1) 



(4.22b) 



-i-w 

where h n is the n b time step. 

For the present equations, one is dealing also with the 
hinge location. Considering the integration of this quantity 
and rewritting the equations (4.22) in a form more suitable 

to iteration, one hasi 



where 



W„ = A ,+ i W h 
n n-l 2 n n 



W = B .+ i W h 
n n-1 2 n n 



h = C ,+ i n h 
n n-i c n n 



A , = W .+ i W , h 
n-1 n-1 2 n-1 n 



B , = W .+ ± W . h 
n-1 n-1 2 n-1 n 



°n-l * %-l + I Vl h n 



(4.23) 

(4.24) 

(4.25) 

(4.26a) 
(4.26b) 
(4. 26c ) 



the values of the acceleration and hinge velocity n n 

are obtained from the equations of motion. 
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For the first mode, the hinge velocity was always zero 
and there is only one equation to be solved, (4.11), which 
turns out to have closed form solution [52]. However, it 
was also integrated numerically to verify the accuracy of 
the numerical calculations and, a perfect agreement was 
achieved. 

For the second mode we have to solve equations (4.15) 

and (4,16). The procedure used in this case was as follows i 

n 

to obtain the values at time step t = E h. begin by deter- 

n t=l 1 

mining the constants (4.26) which will remain unchanged du- 
ring the iterations in the given time step; next assume that 
•• •• • • 

W = W , and h - h and determine (4.23) to (4.25); 

n n-x n n-i 

substitute the calculated values sucessively in (4.16) and 
in (4.15) to obtain and h if these values differ from 

the initially assumed values by differences greater than the 
established convergence criteria, repeat the process until 
satisfactory convergence is attained; when the solution con- 
verged can go to next time step. The solution is considered 
to have converged when it agrees in the fourth decimal place. 
The procedure for the second mode is started by the initial 
conditions* 

»1 = 0 (4.27) 
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1 



(4. 2?b) 



r) 1 = 2- v/2 (4.2 ?c) 

The initial value of acceleration and hinge velocity is 
obtained by substituting (4.2?) in (4.16) and (4.15). 

For the third mode three equations need to be solved, 
for W q , W^, and h , (4.18-20). The general procedure is the 

same as for the second mode except that now at each time 
step we need to achieve convergence in three variables 
instead of only two. 

In Figures 11 and 12 it is shown the time variation 
of all these variables for the first and second modes 
for a value of X = 50. It should be noted that the equa- 
tions solved were in non-dimensional form and the results 
are applicable to the' non-dimensional variables as defined 
before . 

Appendix contains a listing of the computer programs 
used to perform these computations. 

Calculations were made for different values of X so 
as to determine the variation of final deflection with X 
and these results are compared in Figures 8 and 9 with 
the previous solutions and with the experimental results. 
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liI.AUTVO-Pl.AST I G NUMKHI0A1. APPROACH 



The calculations made for the elastic-plastic behaviour 
of the beams, utilized the computer program .JET 3» developed 
by Wu and Witmer [_6^3. 

The program, originally developed to calculate large 
elastic plastic dynamically induced deformations of free and 
restrained, partial and / or complete structural rings, can 
be used to analyse beams since it handles rings of arbitrary 

curvature. 

The program is based on the theoretical development re- 
ported in [6] .It uses the assumed displacement version of 
the finite-element method to model the structure, together 
with the temporal finite-difference approximation. The equa- 
tions of motion for the spatial finite-element formulation 
are derived from the Principle of Virtual Work and D'Alem- 
bert's Principle. The Mises-Hencky yield criteria and it's 
associated flow rule are adopted to describe the elastic 
plastic behaviour of the material. The strain-hardening be- 
haviour is taken into account by using the mechanical sub- 
layer material model. The strain-rate effect is approxima- 
ted by assuming that the uniaxial stress-strain relation is 
affected the strain rate only by a quasi-steady increase in 
the yield stress above the static-test yield stress. The 
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program incorporates the Bernouii i-Euler ( or Kirchoff ) 
hypothesis, excluding transverse shear deformation. 

The behaviour of each finite-clement is characterized 
by four generalized displacements at each node t 



v,w, 



* Fn 



V 

K 



and 




+ 



w 

R 



where v and w are mid-plane displacements in the circun- 
ferential and normal direction, R the radius of curvature 
and H the length cordinate. The displacement behaviour is 
represented by a cubic polynominal in n for both the circun- 
ferential, v , and the normal displacement w. 

Besides elastic restrains JET 3 includes three types 
of nodal displacement conditions namely, symmetry (v=^=0), 
ideally clamped ( v = w = ip = 0 ) and smoothly-hinged 
( v = w = 0 ) . 

As forcing function, the program accepts initial veloci- 
ties and transient externally applied loads. These can be 
described as concentrated at the nodes, uniformly distribu- 
ted or with a sine shape over the elements. 

The equations of motion are solved by applying an appro- 
priate timewise finite-difference operator which provides a 
solution step-by-step in finite-time increments. Two options 
are available for the integration operator t (1) the explicit 
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3-point central difference operator and (:"*) the implicit 
Houbolt operator which ur.es a 4-point backward difference. 

The program is subdivided in four subprograms JET 3A 
to JET 3D which handle four different groups of problems. 

For the present calculations, it was choosen to use 
the Houbolt operator, since it allows the use of larger ti- 
me steps and we are interested in the final deflection of 
the beams analysed. 

However, while for the central difference operator 
there exists a stability criterion for linear dynamic 
systems, which allow the estimation of the appropriate At 
to be used, no such criteria exists for the Houbolt opera- 
tor and the choice of At must be based on numerical ex- 
perimentation. 

Based on the numerical calculations reported in [6] 
for an impulsively loaded beam, it seemed appropriate to 
choose a A t of 5psec. 

While for the central difference method unstable so- 
lution show unreasonably large deflections, for the Houbolt 
operator a gradual degradation of the response may occur 
for large-deflection, non-linear response problems. There- 
fore, initially a calculation was performed also using A t= 
=3psec to assure that a A t of 5|isec was providing the 
converged solution. This was done for specimen 9 of ref. 
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L52] which deforms in the second mode and identical results 
within plotting accuracy were obtained, as shown in Figure 
14. The rest of the calculations were performed using a 
of 5nsec. 

Because of antisymmetry considerations only half of 
the beam was analysed and 10 elements of equal length we- 
re used. The boundary conditions used were fully clamped at 
the support and smoothly hinged at midspan. 

Basically the information the program needs consists 
of the mode locations and the shape of the prescribed ini- 
tial velocity. The initial velocity was given by it's value 
and slope at the nodal stations. As a result, the vicinity 
of the peak velocity was approximated by a higher order 
curve, resulting in higher velocities being introduced. The 
kinetic energy that was given to the program resulted 2 °/« 
higher that in the experiment. 

However, this would correspond to a higher velocity 
and a higher value of X. In Figure 8 the final deflection 
is shown, corrected for new value of X, and Figure 14 
contains the full time history of the response. 

The final deflection is estimated by an average bet- 
ween it's maximum and minimum values after some vibratory 
cycles have elapsed. 
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The same calculations were performed for specimen 10 
of Reference Now 1<> elements of different size wore 

used. I\ smaller one was centered at the peak velocity lo- 
cation and in it’s vicinity finer mesh was used. The same 
type of result was obtained, as can be seen in Figures 8 
and 15 * 
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6. EXPERIMENTAL DETAILS 



The general experimental arrangement used was similar 
to that employed by Jones in several previous studies [28, 
29,49,50,52]. 

The present experiments represent a continuation of 
the studies reported in [,52]. A reader is refered to it for 
a more complete discussion of the experimental details. 

The tests reported here were conducted on fully clamped 

6 2 

beams made of aluminum 6061 T 65 H (E = 10.5x10 lb/in , 

2 4 

P = 0.000251 lb. sec /in ), The beams were nominally 5 in 
long 0,6 in wide and 0.2 in thick. They were loaded with 
layers of sheet explosive 3/8in wide, cut into diamond 
shapes, so as to produce a triangular impulsive velocity 
distribution. 

The impulse imparted to the beam depends not only on 
the weight of explosive but is also influenced by the leader 
arrangement and attennuator. Therefore calibration tests 
have been performed [ 52 ] to determine the necessary cons- 
tants. It resulted that for diamond shaped sheet explosive 
3/8in wide 



W e = w e + 0.03564 gm (6.1) 
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where W is the weight, of the explosive and, for the pro- 
o 

si>nt explosive typo tho specific impulse isi 



I = 0.226.5 lb.sec/gm 

o 



(6.2) 



The total impulse imparted to the beam isi 




(6.3) 



and the peak value of the velocity distribution is 



where d is the length of the beam section which is covered 
with the explosive. 

The ratio between kinetic energy and the maximum amount 
of strain energy which can be absorbed by a structure in a 
wholly elastic manner is estimated byt 



The energy-absorbing characteristics can be measured 
by the dimensionless parameter 



E 



2 

p E V 

o 



(6.5) 



+563 o 

where a Q = 49296 _ 100 r, lb/in . 



Y 




( 6 . 6 ) 



75 



where K q ir, the total initial kinetic energy. 

The permanent transverse displacements were given by 
the difference dial gauge readings before and after the test, 
with the beam still in the clamps. Only the maximum transver- 
se displacements, are shown in Table 1, and are plotted in 
Figures 8 , 9 and 10, together with results from Jones and 

Wierzbicki [52]. 



76 





o\ 


4* 


no 


o 


rH 




i — 1 




On 


O- 




rH 


rH 


U' S 


NO 


O 




CM 




CM 


o 




• 


• 


• 


• 


• 




• 




• 


• 




CM 


CM 


CM 


rH 


CM 




rH 




r ~ 1 


rH 




ON 


vO 


vO 


rH 


CN- 




rH 




on 


UN 


u 


UN 


O- 


O 


iH 


O 




ON 




CM 


ON 


tx3 


♦ 


• 


• 


• 


• 




• 




• 


• 






rH 


i — * 


i — 1 




CM 




0\ 


CM 


• 


(N- 


NO 


O- 


CM 


00 




ON 




cc> 


vO 


W3 


rH 


VO 


CM 


£>- 


UN 




r~^ 




• 


ON 


> r< 


• 


• 


• 


• 


• 




O 




NO 


O 


cd 


UN 


vO 


On 


ON 


h' 




CM 




CM 


CM 














CM -3 


CM 


ON CM ON 


CO CM CO 
















ON r\ON 


00 ONOO. 














ov O 


ON 


vO (S-vO 


o*h a 














rH 


CM 


rH 


CM CM CM 


rH CM rH 


<D 

bQ 

cd X 






















VO 


(N 


ON 


CM 


-3 




On 




O 


On 




CM 


oo 


00 


NO 


i — 1 




UN 




ON 


rH 


0) E 


CM 


CM 


Ht 


rH 


ON 




r— * 




CM 


CM 


> ^ 


• 


• 


• 


• 


♦ 




• 




• 


• 


cd 


































CO 


CO 


ON 


O UN CM 


00 CM ^ 










UN i — 1 


n* O 


ON 


o- 


O 


O-H O 


ON rH -3 










0-3 


0-0 


o 


o- 


ON 


O UN CM 


rH UN O' 


£1 








CM O 


unnC 


o 


o 


vO 


CM CM CM 


0-0-0 


W 

1 

H 








rH CM 


ON CM 


rH 


CM 


rH 


• • • 


iH CM (V 














1 




r 


t t 


f r 
























o 






















CD 






















CO 


on 


O 


rH 


O 


rH 


ON 


C 


ON 


VC O vC 


o' a : a" 


o\ 


o- 


UN 


» — 1 


UN 


H 


O 


00 


O 


On un O 


r — « C r~ \ 


> c 


CM 


-3 


<N 


(X 


CM 


UN 


UN 


UN 


CO On a ' 


UN \ O UN 


•H 


t — 1 


rH 


rH 


rH 


CM 


CM 


CM 


CM 


CM CM CM 


(M CM CM 


























UN 


CM 


UN 


UN 


i — 1 


O 


NO 


Ox 


UN CM UN 


rH O' r— • 


0) E 


00 


-3 


CM 


UN 


ON 


UN 


r-H 


UN 


O UN O 


VO rH vC 


^ B 


ON 




UN 


CM 


0^ 


cm 


CM 


CM 


ON CM ON 


CM C V] CM 




• 


• 


• 


• 


• 


• 


• 


• 


• • • 


# • • 




\C 


-3 


o> 


a: 


UN 




o 




on 


o 


^ — X 


O 


O 


u\ 


ON 


UN 




rH 




o 


UN 


►j c 


O'* 


o- 


un 


CO 


On 




co 




o 


O- 


CM -H 


o* 


o 


Ox 


O* 


O 




Ov 




On 


Ox 




• 


nt 


hJ 




-3 




-3 




-3r 


-3 






ON 


a 


o 


o~ 








CM 


UN 


c 


ON 


ON 


00 


c 


a 




6 




O 


O 


x —i 


Ox 


on 


o 


c 


o 




o 




o 


O 


w 


r-H 

• 


rH 


1 — 1 


CM 


H 




CM 

• 




CM 

• 


CM 




CM 


o 


o 


CM 


o 




NO' 




o 


NO 




Ov 


o 


o 


O 


o 




On 




c 


o 


X C 


ON 


o 


o 


O 


o 




On 




o 


Ox 


*H 


un 


VO 


VO 


vO 


NC 




UN 




VC 






• 


• 


* 


• 


• 




• 




• 


• 


MODE 

No. 


rH 


rH 


1 — 1 


CM 


CM 


ON 


ON 


ON 


E- 

00 * 






















U1 c 


rH 


CM 


on 


-3 


UN 




VC 




o 


cc 



77 



TABLE 1. EXPERIMENTAL RESULTS 



7. DISCUSSION AND CONCLUSIONS 



A study into the higher node dynamic plastic response 
of beams is reported here. Different solutions methods are 
used and results are compared with experiments. 

It is shown that the inclusion of finite-deflections 
improves significantly the infinitesimal deflection solu- 
tions what agrees with previous results ^19-21 , 25,47, 52]. 

The method developed by Jones [_47l account for 
the effect of finite-deflections proved simple to use and 
gave good results. When using the Tresca yield condition, 
as appropriate for metals, reasonable results were obtai- 
ned as can be seen in Figures 8 to 10. The use of a cruder 
description of material behaviour such as the square in- 
teraction curve, without strain-rate effects was much 
simpler to use and rave once more [ 47 , 51 » 52 "] surpri zingly 
pood design estimates o !' final deflection, equations 
( 3. f) , (3.14), (3.18) i as shown in Fi pares H to 10. 

The numerical procedure reported in section 4 which 
satisfies the field equations at every time step is more 
complicated and, for the present problem pave essentially 
the same results as using the previously mentioned method. 
It allowed for the hinge to propagate and it was found 



that the hinge moved only slightly during; the motion and 
the convergence in the numerical iterative process war. very 
rapid what confirms the observe ti ons of Symonds and Chon 
[_4?^]. however, at each time step the moment was computed 
throughout the beam and substantial yield violations were 
found to spread during the motion. 

The el asti.c-plastic numerical analysis is undoutly 
the most compreensive and complex analysis method used, 
however, after being incorporated in a computer code, it 
becomes relatively simple to utilise. As can be seen in 
Figure 9» it gives reasonable predictions, although with 
a tendency of overestimating the deflections. This can 
also be observed in the results reported in where the 

prediction of the final deflection of an impulsively loa- 
ded beam is somewhat higher than the experiment. 

It is apparent from the experimental results that, 
when the initial velocity is arranged with the shape of 
the mode considered, the beam deforms accordingly. Also, 
from experiments reported in (_ 52 ] with uniformly distri- 
buted load, it seem that the beams with different loading 
conditions tend to deform in the mode shapes reported. 
However, it is still to be known how a beam would deform 
when subjected to loading conditions substantially dif- 
ferent from the mode shape, such as a triangular impulses 
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with it's peak not coinciding; with the maximum deflection 
of the mode shape. 

The use of normal modes is a powerfull tool when 
analysing the elastic vibrations of linear systems, the 
only case in which it is exact to use. however, it has 
been used as an approximation in non-linear cases [_65-68], 
It rest to be determined if such concepts could be extended 
to plasticity and the higher dynamic plastic modes used in 
a modal superposition formulation. 
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Figure 1. Free body diarrnm of a beam element :-:hnwW ; 

the positive directions of the indicated 
quantities. For the case of infinitesimal 
deflection;? the axial force N is absent.. 
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Figure 3* First, r->ennd and Third Mode Shaprs 
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Figure 13. Deflected Shape of Specimens 9 and 10 of 
Reference (_ 52j compared with results of 
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EQ A/.(26) REE\52\ 
EQN. 
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APPENDIX 



Contains the results of the analysis described in 
section 4 together with the listing of the computer pro- 
grams used » 

Results for First Mode from Equation (4.11) 

Results for Second Mode from Equations (4.15,16) 
Program Listingi First Mode, Limited Interaction 
Program Listingi Seconde Mode, Limited Interaction 
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FIRST MODE* UNITED INTERACTION! 



DIMENSION B ( 200 ) »BD(200) »BDD(200) ♦ H ( 200 > . T ( 200 ) ♦ A < 200 ) ,BB ( 200 ) 
DIMENSION IT ( 200 ) 

DIMENSION F < 2 1 » 

REAL MM ( 21 ) 

DEFINE CONSTANTS 



I W=5 

RLAM=S0. 

DO 11 1=1 *200 
11 H ( I ) =0 • 0 1 
NREF=5 

DEFINE INITIAL CONDITIONS 



N= 1 

T < 1 ) =0 . 0 
B ( 1 ) =0 • 0 
BD <11=1.0 
BDD(l) =-6./RLAM 

PERFORM THE CALCULATIONS IN TIME STEPS 

WRITE ( IW* 105) 

WRITE ( IW* 1 1 5) RL AM 

10S FORMAT (• «*»*»***«* FIRST MODE * LIMITED INTERACTION ***»*•//) 
115 F ORMAT ( IX t • **«****«« *L AMDA= * »F7.2» •**********•/////) 

1 CONTINUE 
N=N*1 
NN=N-1 

DEFINE THE CONSTANTS FOR EACH TIME STEP 
T (N) =T (NN) 4-H(N) 

A (NN) =BD (NN) ♦BDD(NN) *H(N) /2. 

BB (NN) =B (NN) ♦BD(NN) *H(N> '2. 

IT (N)=0 
B (N) =R (NN) 

INITIATE ITERATION 
I TT = 0 

5 IT (N) =IT (N) ♦ 1 
ITT=ITT*1 
OD=B(N) 

B0D(N)=-6.*(1.*2.*B(N) )/RLAM 
BD (N) =A (NN) ♦BDD(N)*H(N)/2. 

B (N) =BB (NN) ♦BD(N)»H(N) /?_• 

IF (ITT .GT. 5) GO TO 1 
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ERR=DD-B(N) 

IF (ABS(ERR) .GT. 0.0001) GO TO 5 
M = N 

IF ( N .EQ. NREF) GO TO flO 
fl 1 CONTINUE 

IF (N .EQ. 200) GO TO SO 
IF (BO (N) .LT. 0.05) GO TO 45 
GO TO 1 

45 IF (BO (N) .LT. 0.0001 ) GO TO 50 
H(N*1) =0.0025 
GO TO 1 

( 

( CALCULATE the MOMENT AT DISCRETE POINTS ALONG THE BEAM 
( 

80 NREF=NREF*5 

DEFINE THE POINTS TO BE CONSIDERED IN MOMENT CALCULATION 
FR=l./20. 

F (1 ) =0. 

DO 2 I=2»2l 
2 F(I)=F(I-1)*FR 
C CALCULATE THE MOMENT 
DO 40 1=1,21 

MM ( I ) =-RLAM*BDD (N) * (F ( I ) **2. - 1 . ) *F ( I ) /6. ♦ 1 . -2. *F ( T ) 

40 CONTINUE 

WRITE ( IW* 1 1 0 ) N 
WRITE(IW»100) ( F ( I ) , 1 = 1,21) 

WRITE ( I W * 1 0 1 ) ( MM ( I ) , 1 = 1,21) 

100 FORMAT ( 1X.21F6.2) 

101 FORMAT (1X.21F6.2/) 

110 FORMAT ( IX, »*****N=* , 13, »**»**• ) 

GO TO 81 
C 

C MOTION FINISHED. PRINT THE RESULTS 
50 WRITE ( IW*61 ) 

61 FORMAT <*1«) 

WRITE (IW, 105) 

WRITE(IW,115) RLAM 
WRITE ( IW»55) 

55 FORMAT <2X, »TIME» ,RX, *DISPL» ♦ 1 IX, • VEU* , 1 1 X , » ACCEL • , 5X , • IT • ,5X , * STE 
CP* *5X, »DELT* ) 

DO 56 N=1,m 

56 WR I TE ( I W » 60 ) T (N) ,B (N) ,BO (N) ,BDD (N) , I T (N) ,N, H (N) 

60 FORMAT (F6.3,3F15.5,5X, I2.5X, I3,5X,F6.4) 

END 
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SECOND MODE, LIMITED INTERACTION 

DIMENSION R(200> »BD(200> *BDD(200) ,E ( 200 ) ,ED ( 200 ) ,H ( 200 ) 
DIMENSION A (200) ,BB(200) *C ( 200 ) , I T ( 200 ) , I ST ( 200 ) ♦ T ( 200 ) 
DIMENSION r(50) »FF (50) 

HEAL Ml (50) »M2 (50) 

DEFINE CONSTANTS 
I W = 5 

RLAM=50. 

DO 11 1=1*200 
11 H ( I ) =0, 005 
NHEF=5 

DEFINE INITIAL CONDITIONS 



N= 1 

T ( 1 ) =0. 

B ( 1 ) =0 • 

80 ( 1 ) = 1 . 

E(l)=.58S7fl6A376 

BDD(l)=-3.ft(2.-E(N) ) /(R|.AM*E ( 1 ) * ( 1 .-E‘( 1 ) ) ) 
ED ( 1 >=6./<RLAM»E ( 1 ) ) *E ( 1 > tt 8DD< 1 ) 

CALCULATE TIME STEPS 

WRITE ( IW* 105) 

WRITE ( IW« 1 15) RL AM 
1 CONTINUE 
N=N* 1 
NN=N-1 

DEFINE CONSTANTS FOR EACH llME STEP 
T (N) =T (NN) *H (N) 

A <NN> =BD (NN) ♦BOD (NN) »H(N> /2. 

BB(NN)=B(NN) ♦BD(NN) *H(N) /2. 

C(NN)=E(NN) ♦ED(NN)*H(N) /2. 

IT (N)=0 
1ST (N) =0 
B (N) =B (NN) 

E (N) =E (NN) 

BOD (N) =BDD (NN) 

INITIATE ITERATION 
ITT = 0 

5 IT (N) =IT (N) ♦ 1 
ITT= ITT* 1 
DD=B (N) 

BD(N)=A(NN) ♦BDD(N) »H(N) /2. 

B(N)=BB(NN) ♦BD(N) »H(N ) /?• 

IF (ITT ,GT . 5) GO TO 1 
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C BEGIN SUBITEOATI0N 
I NN = 0 

IS 1ST (N) =IST (N) ♦ 1 
I NN= I NN* 1 
D = E <N) 

DB=BDD (N) 

ED(N)=6.*(1 ) / (F <N) *BD <N) *RL AM) *E ( N) *BDD ( N) /BO ( N) 

E (N) =C INN) ♦ED(N) *H(N) /?. 

BOD(N) =-3.» (2.-E (N) (N) ) / ( RLAM*E ( N) * ( l.-E(N) ) ) 

IF (INN .GT.10) GO TO 5 
ERROR = D-E <N> 

IF (ABS(ERROR) .GT. 0.0001) GO TO 15 
ERRO=DB-BDD <N) 

IF (ABS(ERRO) .GT. O.OOOD GO TO 15 
C SUBITERATIONS FINISHED. GO TO NEXT ITERATION 
ERR=DD-B(N) 

IF (ABS(ERR) .GT. O.OOOD GO TO 5 
M=N 

IF (N .EQ. NREF) GO TO 80 
81 CONTINUE 

IF (N .EQ. 200) GO TO SO 
IF (BD(N) .LT. 0.05) GO TO 45 
GO TO 1 

45 IF (BD(N) .LT. 0.0001 ) GO TO 50 
H(N*1) =0.002 
GO TO 1 

80 NREF=NREF*G 

C CALCULATE THE MOMENT AT DISCRETE POINTS ALONG THE BF AM 
C 

C DEFINE THE POINTS TO BE CONSIDERED 
FR=E(M)/20. 

F ( 1 ) =0 . 

DO 2 1=2*21 

2 F ( I ) =F ( 1-1 ) ♦FR 
FF(l)sE(N) 

FFR=(1.-E(M) )/20. 

DO 3 1=2*21 

3 FF(I)sFFd-l) *FFR 
C CALCULATE THE MOMENT 

DO 70 1=1*21 

Ml (I)=(1.-2.*F(I)/E(N) )-RLAM«(BDD(N)-BD(N)*ED(N)/E(N) )*(F(I)** 
C3.0-F ( I ) *E (N) **2.0) / (6.*E (N) ) 

AM1=1 . *E (N) -E (N) *E (N) /?. 

AM2=FF (I ) **3.-3.*FF ( I ) **2. ♦2.*FF ( I ) *4Ml*E (N) »*2.-2.*E (N> 
AM3=BDD(N) ♦BD(N)*ED(N)/(1.-E(N) ) 

H2II)s-d.-FF(I) )/(l.-E(N) ) ♦RLAM*AM3*AM2/(6.*<1.-F<N) ) ) 

70 CONTINUE 

WRITE ( IW* 1 10) N 

WRITE ( IW* 100) (F ( I ) * 1 = 1*21) 

WRITE ( IW* 100) (Ml(l). 1=1*21) 
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WRITE (IW* ion) (FF ( I ) ♦ 1 = 1.21) 

WRITE(IW.IOI) (M2 ( I ) . 1=1.21) 

100 FORMAT ( 1 X . 21F6 • 2 ) 

101 FORMAT ( 1X.21F6.2/) 

105 FORMAT (• ********** SECOND MODE. LIMITED INTERACTION *****•//) 

110 FORMAT ( IX. *****«N=* . 13. ******* ) 

115 FORMAT (IX* »**********i_AMOA=* *F7* 2. •**»*****»*•/////) 

GO TO 81 

C MOTION FINISHED. PRINT THE RESULTS 
50 WRITE ( IW.61 > 

81 FORMAT C • 1 • > 

WRITE ( IW. 105) 

WRITE ( IW. 115) RLAM 
WRITE ( IW.55) 

55 FORMAT (2X. * T ImE • .9X. *HI NGE • ♦ 1 OX t 'DISPL* ♦ llX* • VEL* * 1 1 X ,* ACCEL «. 9X . 
C»HINGEVEL* .5X, * IT* .5X, *Sl • ♦ 5X. • STEP * . 5x . • DEL T* ) 

DO 56 N=1*M 

56 WRITE (IW.60) T(N).E(N),R(N).BD(N) *BD0(N) *ED(N).IT(N).IST(N).N.H(N 
C) 

60 FORMAT (F6,3*5F15.5*7X*l2»5X,I2.5X*lA*5X.F6.4) 

END 
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